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Abstract — A random access scheme for the collision channel 
without feedback is proposed. The scheme is based on erasure 
correcting codes for the recovery of packet segments that are 
lost in collisions, and on successive interference cancellation 
for resolving collisions. The proposed protocol achieves reliable 
communication in the asymptotic setting and attains capacities 
close to 1 [packets/slot]. A capacity bound as a function of 
the overall rate of the scheme is derived, and code distributions 
tightly approaching the bound developed. 



I. Introduction 

Since the introduction of the ALOHA protocol (TJ, several 
random access (RA) schemes have been introduced. Among 
them, some feedback-free RA protocols originally proposed 
in O, J5) re-gained attention in the recent past J4)-||6). In 
|2], the capacity of the so-called collision channel without 
feed-back (CCw/oFB) was analyzed, assuming slot-aligned 
but completely asynchronous users' transmissions. Moreover, 
a simple approach to achieve error-free transmission over 
the CCw/oFB was proposed. The approach of (2) consists 
of assigning different periodic protocol (access) sequences to 
the users. Each sequence defines in which slots each user is 
allowed to access the shared channel. Furthermore, the users 
encode their packets by means of erasure correcting codes. 
The user's packet can be recovered whenever a sufficient 
number of codeword segments are received collision free. 
Hence, by selecting proper protocol sequences, it is possible 
to ensure that a sufficient number of segments per user are 
recovered, even if the beginning of the different protocol 
sequences is unsynchronized. In this way, a capacity equal 
to 1/e [packets/slot] is achieved as M —> oo, where M is the 
number of users accessing the RA channel. The same capacity 
is achieved also in the unslotted case. Although simple, the 
approach of J2) poses some challenges, especially if a large 
(and varying) number of users has to be served 0), H). 

Recently, RA schemes profiting from SIC have been intro- 
duced and analyzed fTl- lfTOl . These schemes share the feature 
of canceling the interference caused by collided packets on 
the slots where they have been transmitted whenever a clean 
(uncollided) copy of them is detected. In J9), iflOl it was shown 
that the successive interference cancellation (SIC) process can 
be well modeled by means of a bipartite graph. The analysis 
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proposed in J9), IflOl resembles density evolution analysis 
of low-density parity-check (LDPC) and doubly-generalized 
LDPC (D-GLDPC) codes over erasure channels fTR-fDl. 
By exploiting design techniques from the LDPC context, a 
remarkably-high throughput (e.g. up to 0.8 [packets/slot]) 
can be achieved in practical implementations. The schemes 
considered in 171- lfTOl assume a feedback from the receiver to 
achieve a zero packet loss rate. Furthermore, the approaches 
proposed in 0, 12) has been recently considered for inclusion 
in the next generation Digital Video Broadcasting Return 
Channel via Satellite (DVB-RCS) standard |[T4l . 

In this paper, we investigate the application of SIC to 
collision channels without feedback. With respect to Q, the 
proposed RA scheme does not require the assignment of 
protocol sequences to users, who select time slots for their 
transmissions in a randomized and uncoordinated fashion. The 
proposed scheme stems from the coded slotted ALOHA (CSA) 
protocol proposed in IflOl . 

Each user splits his generic packet into segments, that are 
encoded through a binary linear (erasure) block code chosen 
on a packet-by-packet basis, according to some probability 
distribution. The resulting encoded segments are then trans- 
mitted in time slots selected with uniform probability out of 
the slots composing the medium access control (MAC) frame. 
If for a packet a sufficient number of encoded segments are 
correctly received, its entire set of encoded segments can be 
recovered by the erasure code. Moreover, the contribution of 
interference of cleaned segments can be canceled from the 
slots in which they caused collisions. It is shown how, iterating 
this process, for a channel traffic below a threshold value, 
arbitrarily low error probabilities can be achieved. We refer 
to this threshold as the capacity of the scheme. By suitably 
selecting the distribution of erasure codes, capacities close to 
1 [packets/slot] can be obtained. It is also shown that there 
is a fundamental trade-off between the rate of the scheme and 
its capacity. An upper bound on the capacity is proved, and 
code distributions tightly approaching it derived. 

II. System Model 

We consider a slotted RA scheme where slots are grouped 
in MAC frames, all with the same length (in slots). Each user 
is frame- and slot-synchronous, and attempts one burst (i.e., 
packet) transmission per MAC frame. Consider M users, each 
attempting the transmission of a burst of time duration T s \ ot 
over a MAC frame of time duration Tf rame . Neglecting guard 
times, the MAC frame is composed of N = Tf rame /T s i ot slots. 
We say that the offered (channel) traffic is G = M/N . 
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Fig. 1 . Model of the access scheme. Each user splits his bursts into k = 2 
segments. User 1 encodes them by a (4, 2) linear block code, users 2 and 
3 by a (3, 2) single parity-check (SPC) code. The darkened squares denote 
parity segments. 



The proposed access scheme works as follows. Prior to 
transmission of a user's burst of time duration T s i ot , the burst 
is divided into k information segments, each of time duration 
T seg = T s \ ot /k. The k information segments are then encoded 
by the user via a packet-oriented linear block code generating 
rih encoded segments, each of time duration T seg . For each 
transmission, the (ri/j, k) code is chosen randomly by the user 
from a set C = . . . , ^n c } of n c candidate codes. For 

h <E {1,2,..., n c } the code % has length n/j, dimension k, 
rate Rh = k/rih, and minimum distance at least 2. Each code 
^ has no idle symbols. At any transmission, each user draws 
his local code from the set C independently of his previous 
choices and without any coordination with the other users. The 
code is drawn according to a probability mass function (p.m.f.) 
A= {Aft,}^f. x which is the same for all users. The MAC frame 
is composed of kN (sub-)slotsQ The coded segments are 
then transmitted by the user over nu slots picked uniformly at 
random. An example for k = 2 is provided in Fig. Q] For the 
special case of k = 1, T sog = T s i ot and each % is a repetition 
code of length rih- The overall rate of the scheme is defined as 
R = k/n, where h := ^h=i ^-h n h is the expected length of 
the code|3 We also define the polynomial A(x) = J2)l=i ^h,x h 
and the polynomial A(a;) = Y^h=i ^h x > where A;,, = 
^■h n h/(z~17=i ^-i n i)- We refer to A(x) as the component code 
distribution. 

We adopt a graph representation of the RA scheme, depicted 
in Fig. H] (related to the same example as Fig. [!)■ Consider 

'The definition of MAC frame as sequence of N slots is instrumental to 
the definition of load G only. The actual minimum unit that can be allocated 
to a segment transmission are the sub-slots. 

2 Note that R does not represent the average code rate R adopted by each 
user. It is easy to prove (Jensen's inequality) that the average rate R = 
y^, h Ajjijjj is lower bounded by R. 
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Fig. 2. Bipartite graph representation of the access scheme of Fig. [T] Subslots 
are represented by SNs, bursts by BNs. Each segment corresoends to an edge. 



M users transmitting over a frame composed of kN slots. 
The situation can be represented by a bipartite graph, Q = 
(B,S,£), consisting of a set B of M burst nodes (one for 
each user), a set S of kN sum nodes (one for each slot), and 
a set £ of edges. An edge connects the i-th burst node (BN) 
to the j-th sum node (SN) if and only if an encoded segment 
associated with the burst of the i-th user is transmitted in the 
j-th slot. The number of edges connected to a BN or SN is 
the node degree. Therefore, a burst encoded via the code ^ is 
represented as a degree-n/j BN, and a slot where d segments 
collide as a degree-e? SN. Moreover, a BN where ^ is used 
during the current transmission is named a BN of type h. 

Each coded segment associated with a type-ft, BN is 
equipped with information about its relevant user and with 
pointers to the other rih — 1 segments of the user0 We assume 
that collisions are always detected at the receiver and that they 
determine a complete loss of information about the colliding 
segments. We also assume that, when a burst is received in 
a clean slot, it is always successfully decoded. Therefore, 
after transmission, we can think of each BN as connected to 
"known" edges and to "unknown" ones, so that some of its 
information segments are known, and the others unknown. At 
a BN of type h, erasure decoding of the code % may recover 
some of the unknown encoded and information segments. It 
is now possible to subtract the interference contribution of the 
newly recovered encoded segments from the signal received 
in the corresponding slots. If d — 1 segments that collided in 
a SN of degree d have been recovered by the corresponding 
BNs, the remaining segment becomes known. The SIC process 
combined with local decoding at the BNs proceeds iteratively, 
i.e., cleaned segments may allow to resolve other collisions. In 
this paper, ideal SIC is assumedQ This procedure is equivalent 
to iterative decoding of a D-GLDPC code over the erasure 
channel, where variable nodes are generic linear block codes 
and check nodes are single-parity-check (SPC) codes. 

3 In practical implementations, the overhead due to the inclusion of pointers 
in the segment header may be reduced by adopting more efficient techniques. 
For fixed k, one may include in the segment header the code index h together 
with a random seed, out of which it is possible to reconstruct (by a pre-defined 
pseudo-random number generator) the positions of the nh segments. 

4 To summarize, our analysis relies on three assumptions: i) Destructive 
(and detectable) collisions; ii) Sufficiently high signal-to-noise ratio (SNR) 
(segments received in clean slots are always successfully decoded); iii) Ideal 
channel estimation (an hypothesis under which, together with the previous 
one, we can assume an ideal SIC). These assumptions simplify the analysis 
without substantially affecting the performance, as shown in (7), (9). 



III. Asymptotic Analysis 

In this section, the evolution of the SIC process for given k 
and G, in the asymptotic case where M — > oo (and N — > oo) 
is overviewed. Most of the results presented in this section 
were developed in [flOl. This serves to define the extrinsic 
information transfer (EXIT) functions for the BN and SN sets, 
used in the proof of our main result presented in Section |IV] 
and to formalize the concept of capacity under iterative SIC. 

We start by briefly recalling the definition of information 
function of a linear block code 03J. Consider an (n, k) linear 
block code c € , where n is the codeword length and k the code 
dimension, and let G be any generator matrix of c € . Then, 
the g-th un-normalized information function of c € , denoted by 
e g , is defined as the summation of the ranks of all possible 
submatrices obtained selecting g columns (with < g < n) 
out of G, regardless their order. 

Next, assume that maximum a-posteriori (MAP) decoding 
is used locally at each BN. At the i-th iteration of the SIC 
process, let Pi-i be the average probability that an edge is 
connected to a SN associated with a slot where a collision 
still persists, before MAP decoding is performed at each BN. 
Moreover, let qi be the average probability that an edge is 
connected to a BN whose contribution of interference on the 
corresponding SN cannot be yet canceled, after MAP decoding 
has been performed at each BN. Then: 

* = i E A * E pt-^ 1 - Pi-ir-'-'iK - mS-t 



h=l 



t=0 



- (t + ljeg^.j =: hipi-i) 



(1) 



Adopting a consolidated terminology in the theory of modern 
error correcting codes, the function /t>(p) is referred to as the 
average EXIT function of the BN set. It is easy to verify that 

hip) = Y2=x X hfl h \p) , where 
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(2) 



is called the average EXIT function of a type-ft, BN, under 
MAP decoding. 

Equation Q] allows to update qi given Pi-\. The dependence 
of pi on qi is instead obtained by observing that the number 
of segments received in a slot follows the distribution 

l / \ M-l 
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where "J; is the probability to receive I segments in the generic 
slot. Defining ty(x) = J2i and letting M — > cxo, yields 



^>(x) = cxp 



(1 - x) 



It follows that pi can be obtained from qi by the relationship 
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Fig. 3. EXIT chart for a regular coded random access scheme employing a 
rate-1/3 repetition code at each BN, characterized by G* = 0.816. 



set. The equations (Q~|) and (O define a discrete-time dynamical 
system qi = with starting point q\ = /b(0), whose 

stability was analyzed in iflOl . Note that the normalized 
offered traffic G is involved in the recursion through (0. 
The asymptotic threshold of the SIC process, denoted by 
G* = G*(C,A), is defined as the sup of the ensemble of 
all G > such that qi — > as i — > oo. In the asymptotic 
setting M -> oo, for all G < G*(C,A) the throughput is 
S = G, i.e., all collisions are resolved. In this sense, G*(C, A) 
represents the capacity of the RA scheme conditional to the 
specific choice of C = {^i, ^i,..., c <a nc } an d 

The recursion defined by (Q]) and Q can be visualized in 
a so-called EXIT chart, which displays fb{p) vs. / s _1 (p). An 
example of EXIT chart for a scheme employing a (3, 1) rep- 
etition code at each BN is provided in Fig. [3] As the iteration 
index i increases, the evolution of the pair of probabilities 
(pi, qi) traces a zig-zag pattern inside the tunnel between the 
two curves. Whenever we operate the scheme below its capac- 
ity, G < G*(C, A), the two curves do not intersect, leaving 
the tunnel open. This lets the probabilities (p, q) achieve the 
(0, 0) point. On the contrary, if the scheme is operated above 
its capacity, G > G*(C,A), the two curves intersect (closing 
the tunnel) in a point (p, q) with p > 0, q > 0, and the iterative 
SIC process converges to a fixed point corresponding to a non- 
zero residual erasure probability. 

IV. Capacity Bound 

In this section, we present our main result along with its 
proof. The proposed proof is based on the EXIT functions 
reviewed in Section [III] and on the Area Theorem, a well- 
known result in the theory of EXIT functions |fl6l . An 
alternative proof (not proposed here for space reasons) can 
be given exploiting algebraic arguments. 



where p(x) = and W(x) = d*(x)/dx. The 

function f s {q) is called the average EXIT function of the SN 



5 Readers with a background in modern coding theory may prefer the 
nomenclature "asymptotic threshold" to "capacity". 



Theorem 4.1: For rational R and < R < 1, let G(R) 
be the unique positive solution of the equation 



V. Capacity- Approaching Schemes 



G = l 



-G/R 



(4) 



in [0,1). Then, the capacity G*(C,A) of the RA scheme 
fulfills 



G*(C,A) < G(R) 



(5) 



for any choice of C = {**?]., ^2, . . . , and A corresponding 
to a rate R. 

Proof: For given C and A, the evolution of the proba- 
bilities (pi,qi) is governed by the recursion pi = f s (qi) and 
Qi = fb(Pi-i) m ® an d IB- f° r all i > 1 and with q\ = /b(0). 
Let's furthermore denote the areas below the BN and the SN 
EXIT functions (over the interval [0, 1]) by 



-4 b 



fb(p)dp and A s = / f s (q)dq 



respectively. A necessary condition for successful decoding is 
to have an open tunnel between the two curves in the EXIT 
chart (see Fig. [3). For all G < G*(C, A), we must have an 
open tunnel in the EXIT chart, which implies 



A, 



As < 1 



(6) 



In particular, © must be satisfied for G — G*(C,A). The 
area below the SN EXIT function (O is given by 



A, 



R a 



R 
G 



(7) 



Moreover, the area below the BN EXIT function (dJ is given by 



A k 



h=l 



X h 



/ b (h) (p)dp=5>,A (8) 



where the second equality holds under the assumption of MAP 
erasure decoding at the burst node, and is due to the Area 
Theorem in 1161 . which states that the area below the MAP 
EXIT function of a binary linear block code without idle 
symbols equals its code rate. By incorporating (0 and ([H) 
in ^ we obtain 
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(9) 



Observing that X h = n h A h / J2i( n i^i) = n h k h /n, we have 
"Y^ukXh/rih = k'}2 h h^i l /n = k/n = R. Thus, <j9j reduces to 
R + {R/G)e~ G / R < R/G which, for G = G*(C, A) yields 



G*(C,A) < 1 



G*(C,A) 



(10) 



Defining G(R) as the unique solution in (0, 1] of ©, we 
have R(G) = -G/log(l - G). Moreover, from ([Toll we 
have R(G*) > -G*/log(l - G*). Since the function y = 
— m/log(l — x), x G [0, 1), is monotonically decreasing, we 
obtain (0. □ 

Note that, while G*(C, A) depends on R through C and A, 
its upper bound G(R) depends solely on R. 



Theorem 14. 1 1 establishes an upper bound on the capacity of 
the proposed RA scheme. In this section, we show that this 
bound can be indeed approached. To do so, for a given k, 
a given set C = \f€\ , ^2, - - . , ^ } of component codes, and 
a given target rate R, we generated by differential evolution 
(DE) lUTI optimization the distribution A which maximizes 
G*(C, A). In order to limit the search space, we focused on 
schemes based on codes of moderate-low length. 

For example, for the case of k = 1 and R = 0.2, we limited 
nh to 30 (i.e., repetition codes with rate down to 1/30 have 
been considered). In this case, we obtained the distribution 

Ai(x) = 0.494155x 2 + 0.159085x 3 + 0.107372x 4 
+ 0.070336x 5 + 0.045493a; 6 + 0.019898a; 7 + 0.024098a; 11 
+ 0.008636a; 12 + 0.005940a; 13 + 0.008749a; 15 + 0.002225a; 18 
+ 0.001261a; 20 + 0.002607a; 22 + 0.008092a; 23 + 0.00 2 2 8 7a; 24 
+ 0.012274a- 25 + 0.00 2 5 30a; 26 + 0.003094a; 27 + 0.002558a; 28 
+ 0.005891a; 29 + 0.0 1 34 1 9a; 30 , 

where Ai(x) = Ylh ^hX h and % is the (h, 1) repetition code. 
The associated threshold is G* = 0.977. The corresponding 
point on the G x R plane is reported in Fig. |4]and compared 
with the bound given by Theorem 1. On the same plot, 
other distributions, denoted by A* (a;) for i <E {2, . . . , 5}, are 
reported for k = 1. Whereas for low rates R repetition-based 
scheme approach tightly the bound, for rates close to 1/2 they 
show visible losses. The distribution As(x) = 0.8a; 2 + 0.2a; 3 
(obtained by allowing a maximum nu of 5) presents a rate 
R = 5/11 ~ 0.45 and attains a threshold G* = 0.625, whereas 
G(5/ll) ~ 0.843. This effect if somehow expected, since 
in the limiting case of R = 1/2 each user employs a (2, 1 
repetition code, and the corresponding threshold is G* = 0.5 

As the rate R grows, it is convenient to adopt codes 
with k > 1. To this purpose, we designed further schemes 
where each ^ € C is maximum distance separable (MDS)Q 
It is assumed that each burst node adopts a sub-optimum 
(in the symbol-wise MAP sense) decoding approach, which 
consists of decoding its lost segments only if the number of 
its collision-free received segments is at least k. Under this 
assumption, the EXIT function of a BN using an (rih, k) MDS 
codes is given by 



fi h \p) = Y.( nh l l \ 

1=0 ^ 



-l-l 



(11) 



We define in this case the distribution A(x) = J^h ^-h% \ 
where % is an (h + k, k) MDS code. 

For k = 2 we designed a distribution Kq(x) = 0.276023x4- 
0.366641a; 2 + 0.127979a; 3 + 0.229357a; 7 characterized by 
R = 0.4 and G* = 0.83. For the same rate, the best found 
repetition-based (k — 1) distribution for a maximum rih set 
to 10, denoted by A 4 (a;), achieves G* = 0.79. Moving to the 

6 For k = 1, R = 1/2 represents a limiting case because, for rates R > 
1/2, the rate-1 repetition code (i.e., no coding) must be included in the set 
C. In this case, it is easy to prove that the threshold G* would drop to 0. 

'imposing limits on n^, this assumption is realistic. For instance, (gener- 
alized) Reed-Solomon codes on finite fields of moderate order may be used. 
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Fig. 4. Upper bound to the capacity vs. rate R. Thresholds G* are reported 
for selected distributions. Distributions Ai(x), i = 1, . . . , 5 (*) are based on 
repetition codes (k = 1). Distribution Aq(x) (□) is based on MDS codes 
(k = 2). Distributions A{(x),i = 7,8,9 (A) are based on MDS codes 
(k = 3). Distributions Ai(x),i = 10,11 (~) are based on MDS codes 
(k = 4). Distributions based on (k + 1, k) SPC codes are also displayed (o). 



moderate-rate regime, we heuristic ally found that the bound 
can be better approached by adopting distributions based on 
higher code dimensions (see Fig. |4). We conjecture that the 
bound can be approached even more tightly for larger code 
dimensions (fc). In Fig. [4] the thresholds achieved by regular 
schemes based on SPC codes of increasing rates are provided. 
In has been proved in iflOl that for these schemes the threshold 
admits the close form G* = (k + 1) — 1 . As k grows, the 
rate of these scheme, R = k/(k + 1) approaches 1 and the 
corresponding threshold tends to 0. For large k, the scheme 
tends to operate close to the capacity bound for high rates. 

Fig. shows the packet loss rate (PLR) achieved by the 
scheme employing the distribution Ai(cc). The results have 
been derived via Monte Carlo simulations for MAC frames of 
size N — 5000, 1000 and 500 slots, and are compared with 
the capacity of the scheme, G* = 0.977. For the N = 5000 
case, a PLR close to 2 • 10~ 3 is achieved at a channel traffic 
G = 0.94 [packets/slot], only 0.05 [packets/slot] away from 
the bound established by Theorem|4j](~ 0.99 [packets/slot]). 



VI. Conclusions 

We introduced a RA protocol for the CCw/oFB which 
achieves large efficiencies with mild coordination demands. 
The scheme is based on the use to erasure correcting codes to 
recover packet segments that are lost in collisions, and on SIC 
to resolve collisions. The proposed protocol achieves reliable 
communication in the asymptotic setting and attains capacities 
close to 1 [packet /slot]. A simple capacity bound as a function 
of the code rates employed to encode the segments has been 
proved. The derived bound can be approached by means of 
judiciously designed component code distributions. 



Fig. 5. PLR for the RA scheme based on the distribution Ai(x). N = 
5000, 1000, 500, maximum iteration count set to 100. 
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